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Abstract 

Dynamical simulations of polaron transport in conjugated polymers in the presence of an exter- 
nal time-dependent electric field have been performed within a combined extended Hubbard model 
(EHM) and Su-Schrieffer-Heeger (SSH) model. Nearly all relevant electron-phonon and electron- 
electron interactions are fully taken into account by solving the time-dependent Schrodinger equa- 
tion for the vr-electrons and the Newton's equation of motion for the backbone monomer displace- 
ments by virtue of the combination of the adaptive time-dependent density matrix renormalization 
group (TDDMRG) and classical molecular dynamics (MD). We find that after a smooth turn-on of 
the external electric field the polaron is accelerated at first and then moves with a nearly constant 
velocity as one entity consisting of both the charge and the lattice deformation. An ohmic region 
(3 mV/A < Eq < Q mV/A) where the stationary velocity increases linearly with the electric field 
strength is observed for the case of C/=2.0 eV and V=1.Q eV. The maximal velocity is well above 
the speed of sound. Below 3 mV/A the polaron velocity increases nonlinearly and in high electric 
fields with strength Eq > 10.0 mV/A the polaron will become unstable and dissociate. The rela- 
tionship between electron-electron interaction strengths and polaron transport is also studied in 
detail. We find that the the on-site Coulomb interactions U will suppress the polaron transport 
and small nearest-neighbor interactions V values are also not beneficial to the polaronic motion 
while large V values favor the polaron transport. 
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I. INTRODUCTION 



Charge transport properties in conjugated polymers have attracted sustained attention 
from both academic and industrial researchers, since it was discovered in the 1970s that 
the electrical conductivity of trans-polyacetylene (PA) can be improved significantly by 
doping with strong electron acceptors or donors^, 3,j3|. Many electronic devices have been 




fabricated based on these conjugated polymers. 

All the conducting polymers can generally be sorted into two classes. The first class 
is trans-polyacetylene, in which there exists a twofold degeneracy of ground state energy 
distinguished by the positions of the double and single bonds. The degenerate ground state of 
trans-polyacetylene leads to solitons 8|] as the important excitations and the dominant charge 
storage species, which take the form of domain walls separating different districts of opposite 
single and double bonds alternation patterns. The second class of conducting polymers 
is given by all the other systems except trans-polyacetylene, in which the ground state 
degeneracy is weakly lifted. So polarons and confined soliton pairs (bipolarons) 8|] instead 
of solitons are the important excitations and the dominant charge storage configurations 
for this class of conducting polymers. Nowadays it has been widely accepted that these 
quasiparticles (solitons, polarons and bipolarons) are the fundamental charge carriers in 
conducting polymers. The studies of the dynamics of these charge carriers are therefore of 
great interest for the purpose of modulating or devising new organic electronic materials 
based on conjugated polymers. 

So far, there have been many extensive theoretical studies on polaron dynamics in conju- 



gated polymers 



(SSH) model 
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Q, S] based on the Su-Schrieffer-Heeger 



|, an improved Hiickel molecular orbital model, in which vr-electrons 
are coupled to distortions in the polymer backbone by the electron-phonon interaction. 
Rakhmanova and Conwell [o], [l^ have considered the motion of a polaron under low and 
high electric fields respectively by using an adiabatic simulation in which the electronic en- 
ergy is treated within the Born-Oppenheimer approximation. Their results show that the 
polaron can form and move in its entity in low electric fields but dissociates in high fields. 
The highest electric field strength under which the polaron can sustain is about 6 mV/A. 
Basko and Conwell ll| have also analyzed the dynamics of a polaron and found that the 
speed of a steadily moving polaron cannot exceed the sound speed (Vs). These conclusions 
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are consistent with those from the numerical results presented by Arikabe et al 121 but 

n n 

inconsistent with those given by Johansson and Stafstrom [13|, [1^, who performed a nona- 
diabatic simulation in which transitions between instantaneous eigenstates are allowed. In 
Johansson and Stafstrom's work, the maximum velocity (about ^Vs) is reached at a high 
electric field strength (~ 3.5 mV/A) and at even higher electric field strengths the polaron 
will become unstable and dissociate. 

In the SSH model only the electron-lattice interactions are considered and the electron- 
electron interactions are totally ignored. In order to improve this situation, Di et al have 
considered extended Hubbard model (EHM) combined with the SSH model and investigated 
the dynamics of polarons within this model at the level of an unrestricted Hartree-Fock 
(UHF) approximation. 23|] They found that the localization of the polaron is enhanced and 
the stationary velocity of the polaron is decreased by both the on-site repulsions U and the 
nearest-neighbor interactions V . Furthermore, they found that the local extremum of the 
stationary velocity of the polaron appears at ?7 ~ 2V . However, in recent years a lot of theo- 
retical calculations of the static properties of conjugated polymers have shown that the elec- 



tron correlation effect plays a very important role in determining the behavior o 



carriers in conjugated polymers. 
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For example, electron correlation is essential in obtaining correct pictures of the spin polar- 



izations in trans-polyacetylene. 
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Therefore, theoretical simulations for polaron dynamics 
with electron correlations considered, which step beyond the SCF level, are highly desirable. 
However, the calculations for large polymer systems with traditional advanced electron- 
correlation methods such as configuration interaction method (CI), multi-configuration self- 
consistent field method (MCSCF), many-body MoUer-Plesset perturbation theory (MPn) 
and coupled cluster method (CC) are currently still not feasible due to the huge compu- 
tational costs. Fortunately, the ada pti ve time-dependent density-matrix renormalization 



group (TDDMRG) method 



41 



42, 



431] can be used instead. In the context of ID cor- 



related electronic and bosonic systems, the adaptive TDDMRG has been found to be a 
highly reliable real-time simulation method at economic computational cost, for example in 



the context of magnetization dynamics 



44| , of spin-charge separation 45|, |46| , or far-from 



equilibrium dynamics of ultracold bosonic atoms 47|. Recently, Zhao et al have performed 
adaptive TDDMRG simulations for polaronic transport within a combined model of SSH 
and Hubbard model to take the on-site Coulomb interactions into account and found that 
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that the velocity of the polaron is suppressed by the on-site Coulomb interaction U. 48 1 

In this paper, combining the SSH and extended Hubbard model to take both the on- 
site Coulomb interactions and nearest-neighbor electron-electron interactions into account, 
we simulate the motion of a polaron in conjugated polymers under an applied external 
electric field by using the adaptive TDDMRG for the vr-electron part and classical molecular 
dynamics (MD) for the lattice backbone part. The aim of this paper is to give an exhaustive 
picture of polaron transport in conducting polymers at a theoretical level with all relevant 
electron-electron interactions and correlations included and to show how the on-site Coulomb 
interactions U and nearest-neighbor electron-electron interactions V influence the behavior 
of polaron transport in conducting polymers. 



II. MODEL AND METHODOLOGY 



We use the well-known and widely used SSH Hamiltonian 2l|, |22| combined with the 
extended Hubbard model (EHM), and include the external electric field by an additional 
term: 



H{t) = Hei + HE(t) + Hi 



llatt 



(1) 



This Hamiltonian is time-dependent, because the electric field E[t) is explicitly time- 
dependent. 

The TT-electron part includes both the electron-phonon and the electron-electron interac- 
tions, 



Hel — — ^ tn,n+l{c^+l^a'^n,a + h.C.) 



(2) 



2''v''n+l,o-'''"+l.o-' 2' 



where tn,n+i is the hopping integral between the n-th site and the {n + l)-th site, while 
U is the on-site Coulomb interaction and V denotes the nearest-neighbor electron-electron 
interaction. Because the distortions of the lattice backbone are always within a certain 



limited extent, one can adopt a linear relationship be 
lattice displacements as tn,n+i = ~ «(wn+i — Un) [21 



;ween the hopping integral and the 



22], where to is the hopping integral 
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for zero displacement, Un the lattice displacement of the nth site, and a is the electron- 
phonon coupling. 

Because the atoms move much slower than the electrons, we treat the lattice backbone 
classically with the Hamiltonian 

n n 

where K is the elastic constant and M is the mass of a site, such as that of a CH monomer 
for trans-polyacetylene. 

The electric field E{t) directed along the backbone chain is uniform over the entire system. 
The field which is constant after a smooth turn-on is chosen to be 

fEoexp[-(t-Tc)VT^], fort<Tc, 
E{t) = I (4) 

[Eo, ioit>Tc, 

where Tq, Ty/ and Eq are the center, width and strength of the half Gaussian pulse. This 
field gives the following contribution to the Hamiltonian: 

HE{t) = \e\ ^{na + Un){c^^„Cn,a - \)E{t) (5) 

71, cr 

where e is the electron charge and a is the unit distance constant of the lattice. The model 
parameters are those generally chosen for polyacetylene: to=2.5 eV, a=4.1 eV/A, K=21 
eV/A^, M=1349.14 eVfsVA^ a=1.22 A. [22] The resuhs are expected to be qualitatively 
valid for the other conjugated polymers. The values of Tc and Tw (Tc=30 fs and Tw=25 
is) are taken from the paper of Fu et al. [3] 

For the purpose of performing real-time simulation of both the evolution of quantum vr- 
electron part and the classical movement of the chain backbone, we adopt a newly developed 
real-time simulation method in which classical molecular dynamics is combined with the 
adaptive TDDMRG. The main idea of this method is to evolve the vr-electron part by the 
adaptive TDDMRG and move the backbone part by classical MD iteratively. Details about 
this method can be found in recent papers Q, |49 |. 



III. RESULTS AND DISCUSSION 



To investigate the dynamic properties of polarons, we simulate the polaron transport 
process in a single model conjugated polymer chain under a uniform external electric field. 
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FIG. 1: Time evolution of net charge distribution in a polymer chain with the transport of a 
polaron under an external electric field. (The polaron moves from the left to the right.) (£"0=3.0 
mV/A, U=2.0 eV, V=1.0 eV) 



A polaron is different from a charged soliton in that it involves an unpaired electron and 
thus a nonzero 5 = 1/2. In our calculations, we simulate a model chain containing = 100 
monomers and = 99 vr-electrons to present a polaron defect without degenerate ground 
states. The polaron is initially centered around site 30 and site 31 through imposing a 
constraint of reflection symmetry around the center between site 30 and site 31. Then, the 
dynamics of polaron transport in 160 femtoseconds (fs) is simulated by virtue of classical 
MD combined with the adaptive TDDMRG. 



A. General picture of polaron transport 

Firstly let us show the general time evolution picture of polaron transport in a conducting 
polymer. The time evolution of charge density and the staggered bond order parameter 
r„ = (— 1)"(2m„ — Un+i — Mn-i)/4 in a polymer chain with a polaron defect are shown in 
Fig. [1] and Fig. [21 It can be clearly seen that both the charge density and the geometrical 
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FIG. 2: Time evolution of the staggered bond order parameter r„ in a polymer chain with the 
transport of a polaron under an external electric field. (The polaron moves from the left to the 
right.) (£;o=3.0 mV/A, U=2.0 eV, V=1.0 eV) 



distortions don't localize at one monomer. On the contrary, the polaron defect spreads over 
a delocalized region with a length of tens of sites. One can also clearly see that during the 
entire time evolution process the geometrical distortion curve and net charge distribution 
shape for the polaron defect always stay coupled and show no obvious dispersion. This 
implies that the polaron defect is an inherent feature and the fundamental charge carrier in 
conducting polymer. In Fig. [21 we can also see that a long lasting oscillatory "tail" appears 
behind the polaron defect center. This "tail" is generated by the inertia of those monomers 
to fulfill energy and momentum conservation; they absorb the additional energy, preventing 
the further increase of the polaron velocity after a stationary value is reached. 

In order to evaluate the velocity of polaron transport, it is necessary to know the center 
postions of a polaron at different times. The center position of a defect can be derived from 
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FIG. 3: Temporal evolution of the center position of a polaron under an external electric field 
derived from both the charge density and the bond order parameter. (£'o=3.0 mV/A, U=2.0 eV, 
V=1.0 eV) 

the charge density distribution picture or the bond length alternation pattern as 



N{7r + e,^g{t))/27T, 



if cos6'c,g > and sin6'c,c, > 0; 
if cosOc^g < 0; 



(6) 



N{2tt + 6'c,g(t) )/27r, otherwise, 

where 6c,g is defined according to the charge density pn or the staggered bond order parameter 
r„ as 

E„PnWsin(27m/iV) 

(7) 



arctan 



arctan 



E„Pn(t)cos(27m/iV)' 
En(^n,(t)-r°)sin(27™/iV) 



E„(^n(t)-rO)cos(27m/iV)' 

in which r° is the dimerized lattice displacement of the pristine chain without defects. 

From Fig. [Hand Fig. [21 one can also find that, after the external electric field is smoothly 
turned on, the polaron is accelerated at first and then moves with a nearly constant velocity 
as one entity consisting of both the charge and the lattice deformation. We calculate the 
postions of the center of a polaron during the transport process according to Eq. [6]and show 
the temporal evolution of the polaron center position in Fig. [3l It can be clearly seen that 
the center position of the polaron increases nearly linearly with the increase of time after the 
initial short-time acceleration. This implies the polaron defect is transported with a nearly 
constant velocity and supports our observation from Fig. [T] and Fig. [21 Therefore, only 
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FIG. 4: The stationary velocity of a polaron Vp as a function of the electric field strength. (U=2.0 
eV, V=1.0 eV) 

the stationary velocity will be considered for polaron transport in the following discussions. 
Meanwhile, one also finds that the lattice distortion center always stays well together with the 
charge density center. But it should also be noticed that, under high electric field strength 
the deviation of the lattice distortion center position from the charge density center position 
may become much larger because the charge density moves much faster under high field 
and the lattice distortion cannot catch up with the charge density. This separation between 
the charge center and lattice distortion center may lead to the dissociation of a polaron. 
The dissociation of a polaron under high electric field strength will be discussed in the next 
section. 
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B. Relationship between polaron transport velocity and electric field strength 

The dependence of the stationary velocity of a polaron Vp on the electric field strength is 
shown in Fig. HI We find that Vp increases with increasing electric field strength. Apparently, 
the polaron velocity can easily exceed the sound speed {Vs = \/4:K / Ma/2=Q.\b A/fs) and 
the maximum velocity (about SV^) is reached at a high electric field strength (~ 9 mV/A). 
Moreover, a very interesting behavior of function of the electric field strength is 

observed. An ohmic region where Vp increases nearly linearly with the electric field strength 
is found. This ohmic region extends approximately from 3 mV/A to 9 mV/A for the case 
of t/=2.0 eV and V=1.Q eV. Below 3 mV/A the polaron velocity increases nonlinearly and 
when Eq > 10 mV/A the polaron becomes unstable and the charge will be decoupled from 
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FIG. 5: Temporal evolution of the center position of a polaron under different external electric 
field strengths. (U=2.0 eV, V=1.0 eV) 

the lattice distortion. 

We illustrate the temporal evolutions of the polaron center postion under different exter- 
nal electric field strengths. We can find that the lattice distortion center stays well together 
with the charge center under low electric field strengths. However, the deviation of the 
lattice distortion center from the charge density center becomes much larger under higher 
electric field strengths because the charge density moves much faster under high field and 
the lattice distortion cannot catch up with the charge density. Under the electric field of 
10 mV/A, the dist ance between the two centers becomes even larger than several unit dis- 
tances, implying that the polaron has been completely dissociated. This dissociation can 
be clearly seen from the time evolution pictures of charge density and the staggered bond 
order parameter in a polymer chain with a polaron defect under a external electric field of 
10 mV/A in Fig. [6] and Fig. [71 While the charge density accelerates to a high speed as shown 
in Fig. El the original lattice distortion stays far behind the charge density with a distance 
of more than tens of unit distances. The charge density will induce new lattice distortion 
in the nearby region around the charge density center, while the original lattice distortion 
can hardly move further, as illustrated in Fig. [71 Therefore, the polaron defect has become 
completely dissociated under the high electric field. Apparently, our calculated threshold 
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FIG. 6: Time evolution of net charge distribution in a polymer chain with the transport of a 
polaron under a high external electric field. (The polaron moves from the left to the right.) 
(£;o=10.0 mV/A, U=2.0 eV, V=1.0 eV) 

field strength of around 10 mV/A for polaron dissociation is in agreement with the value of 
9.5 mV/A obtained by SSH-EHM/UHF calculations by Di et al 23|] and much closer to the 
experimental determined value of about 15 mV/A 50] compared to the value of 6 mV/A 



calculated 
Stafstrom 



Rakhmanova and Conwell 0, [l^ and 3.5 mV/A calculated by Johansson and 
13|,[l4|. The large deviations of the calculated threshold field strength from the 
experimental value by the latter two groups are due to the fact that they considered only 
the SSH model and considered no electron-electron interactions in their studies, and the 
agreement of SSH-EHM calculations with experimental results shows again that the inclu- 
sion of electron-electron interactions is vital for the theoretical studies of charge carriers in 
conducting polymers. Another thing that should be noticed is that our U value (2.0 eV) 
and V value (l.OeV) are not the standard values for any special kind of polymers. Here 
we use these two values only for the purpose to show how a polaronic transport will gen- 
erally evolve and whether the inclusion of electron-electron interactions is important for a 
reasonable theoretical simulation of the polaronic transport process. The detailed discussion 
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FIG. 7: Time evolution of the staggered bond order parameter r„ in a polymer chain with the 
transport of a polaron under a high external electric field. (The polaron moves from the left to the 
right.) (£;o=10.0 mV/A, U=2.0 eV, V=1.0 eV) 

of how the on-site repulsions U and the nearest-neighbor interactions V will influence the 
polaronic transport will be shown in the next sections. 

C. Influence of on-site repulsions U on polaron transport velocity 

In order to study the influence of electron-electron interactions on polaron transport, we 
focus on the stationary velocity of a polaron Vp calculated with different electron-electron 
interaction strengths under the constraint that the other parameters are fixed. As should 
be noticed that, the real conjugated polymer is with weak interactions. Strong interactions 
will lead to too strong charge polarizations which are unrealistic. Considering that we are 
only focusing on the study of real conjugated polymer system, in this work we adopt only 
the weak-interaction parameters {U < 2t,V < t). 

Firstly, we study the Hubbard model with only on-site Coulomb interactions f/, i.e., 
\^ = in the EHM. The dependence of the polaron stationary velocity Vp on the on-site 
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FIG. 8: The stationary transport velocity of a polaron as a function of different U values. {Eq=3.0 
mV/A, V=0.0 eV) 
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FIG. 9: The staggered bond order parameter of a static polaron for several different U values. 
{V=0.0 eV) 

Coulomb interactions U is displayed in Fig. [HI We find that Vp decreases monotonically with 
the increasing U value and its maximum value is achieved at U=0.0 eV. This can be easily 
understood because the variation of the charge carrier transport velocity is strongly related 
to the delocalization level of the charge carrier defect. The lattice tends to be occupied 
by one electron per site when the on-site Coulomb repulsion U increases. Therefore an 
increasing U will lead to more localized charge density and a smaller defect width of the 



13 



polaron and accordingly restrain the polaron transport. This analysis is supported by our 
calculated charge density width values of a static polaron W^- Wc is calculated according to 
the following formula: 

W, = {Y,[n-X,fpnf'\ (8) 

n 

Through calculations we find the polaron width Wc also decreases monotonically with an 
increasing U value, from 10.10 (t/=0.0 eV) to 9.68 (f/=1.5 eV) and then to 7.62 (f/=4.0 eV). 
This sequence is in agreement with that of the polaron stationary velocity Vp and supports 
our analysis of the relationship between the polaron transport velocity and the delocalization 
level of a polaron. Actually, the change of the polaron delocalization can also be directly 
viewed through the geometrical picture of the static polaron defect. In Fig. [HI the staggered 
bond order parameter r„ of a static polaron calculated by different U values is shown. It 
can be clearly seen that the lattice dimerization is enhanced and the polaron width becomes 
narrower while the U value increases. All these pictures show that the polaron transport is 
suppressed by the on-site Coulomb repulsions U . The monotonic decrease picture of polaron 
transport velocity with on-site Coulomb repulsions U observed by us is in agreement with 



Di et a/'s UHF calculations |23| and Zhao et aUs recent adaptive TDDMRG results 48|], but 
not in accordance to our recent studies on charged soliton transport, in which we found that 
charged soliton transport velocity is non- monotonic in U 49|]. The different behaviors of 
charge carrier transport with U increasing is due to the different characteristics of charged 
solitons and polarons. It was found that the polaron defect is more delocalized than the 
charged soliton defect and the height of the charge density peaks in polarons is only roughly 
one half that of the charge density peaks in charged solitons. 39| In a more localized case 
with large charge densities, as in a charged soliton, small U may favor the charge carrier 
transport because the increasing on-site repulsion for large charge densities can make the 
electron (or hole) to hop more easily to the next site. Therefore the different behaviors of 
charge carrier transport with U increasing for charged solitons and polarons are reasonable. 



D. Influence of nearest-neighbor interactions V on polaron transport velocity 

Secondly, we also study the influence of nearest neighbor electron-electron interactions V 
on the polaron transport process. The dependence of the stationary velocity of a polaron Vp 
on the V values is displayed in Fig. [TUl where U values are supposed to be frozen to be zero. 
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FIG. 10: The stationary transport velocity of a polaron as a function of different V values. {Eq=3.0 
mV/A, U=0.0 eV) 




FIG. 11: The charge density of a static polaron as a function of the site index calculated with 
different V values. (U=0.0 eV) 

This assumption is not realistic because the on-site Coulomb interactions U are normally 
much stronger than the neighbor electron-electron interactions V; we make this assumption 
only for the purpose of studying of the influence of V on the polaronic transport process 
without the effect of U. As can be seen in Fig. [TOl the situation is completely different 
from the case of an on-site Coulomb interaction. Interestingly, Vp is non-monotonic in V: it 
decreases to a shallow minimum at V ^ 1.2 eV gradually and then increases again. In fact, 
similar to the the case discussed above with varying U, the variation of Vp with V is also 
strongly related to changing the delocalization of the polaron defect. For small V, larger 
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FIG. 12: The staggered bond order parameter r„ of a static polaron for several different V values. 
(C/=0.0 eV) 

positive charge densities will be induced in the central part of a polaron defect while the 
induced negative charge densities at nearest neighbor sites are still relatively small, as shown 
in Fig. [m Therefore the polaron defect tends to be more localized and consequently the 
polaron transport becomes more difficult. So, Vp decreases while the value of V increases 
from 0.0 eV to 1.2 eV as displayed in Fig. [101 However, when nearest neighbor electron- 
electron interactions V increase further and become dominant, very large charge polarization 
will be induced. Therefore, the electron-hole attraction between opposite charge densities 
at nearest-neighbor sites will contribute much more significantly to the polaron system and 
favor the hoppings of the accumulated electrons (or holes) in the central part of a polaron 
defect to the neighbor sites. As shown in in Fig. [H], the two positive charge density peaks 
are separated further when V increases from 1.2 eV to 1.4 eV. This change leads to a 
more delocalized polaron defect and accordingly the increase of Vp. Our calculated polaron 
charge density width results, namely that Wc decreases gradually from 10.10 (V^=0.0 eV) 
to 9.37 {V=1.2 eV) and then increase gradually to 9.89 (1^=1.4 eV), support our analysis 
well. In order to directly view the change of delocalization level of the polaron defect through 
geometrical pictures, we also show the staggered bond order parameter r„ of a static polaron 
calculated with different V values in Fig. [121 It is clearly shown that the lattice dimerization 
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FIG. 13: The stationary transport velocity of a polaron as a function of U for fixed values of V. 
{Eq = 3.0 mV/A) 



is enhanced and the polaron width becomes narrower with V for small V . We also find that 
the polaron becomes more delocalized with a smaller r„ minimum while V increases further 
from 1.2 eV to 1.4 eV. The change of polaron delocalization illustrated from the geometrical 
picture is in accordance to that has been shown in the charge density picture, verifying 
that small V suppresses the polaronic transport while large V favors the polaron transport. 
This non-monotonic picture of polaron transport velocity in V is similar to that has been 



observed in charged soliton transport 



49| 



E. Influence of both the on-site Coulomb interactions U and the nearest-neighbor 
interactions V on polaron transport velocity 

Furthermore, we consider the realistic case of conducting polymers in which both the on- 
site Coulomb interactions U and the nearest-neighbor interactions V are taken into account. 
Fig. [13] shows the stationary transport velocity of a polaron Vp as a function of U for fixed 
values of = 0.2 eV, V = 0.5 eV, ^=1.0 eV and V = 1.5 eV. Generally, the Vp behavior 
with variated U and fixed V values is quite similar to that has been illustrated in Fig. [S] in 
which V is frozen to be 0.0 eV. It can be found that Vp decreases monotonically with U. 



17 



Meanwhile, we can also find that, Vp seems to decrease more rapidly in large V cases than in 
small V cases for small U . Larger nearest-neighbor electron-electron interactions V will lead 
to larger charge polarizations and of course the increasing on-site Coulomb repulsions U will 
suppress the transport of a polaron associated with larger charge polarizations more signif- 
icantly comparing to the case with smaller charge polarizations. Therefore the suppression 
of polaron transport by the on-site Coulomb repulsions U is more remarkable in large V 
cases than in small V cases. The monotonic decrease picture of polaron transport velocity 
with on-site Coulomb repulsions U observed by us is opposed to the non-monotonic picture 
obtained by Di et al through UHF calculations 23|. They found a local extremum of the 
stationary velocity of the polaron can be achieved at f/ ^ 2V . The difference between the 
results of the adaptive TDDMRG and UHF calculations is due to the fact that the latter 
method ignores the important electron correlation effect. Previous theoretical static studies 
have found that neglecting electron correlation effects will lead to an underestimation of 



34 



35 



31 In a 



the defect width as well as an overestimation of charge polarization, 
more localized case with large charge densities, presented by UHF calculations without the 
inclusion of electron correlation effect, introducing a small U may favor the charge carrier 
transport because the increasing on-site repulsion for large charge densities can make the 
electron (or hole) to hop more easily to the next site. Therefore, the difference between TD- 
DMRG calculations and UHF calculations shows again that the electron correlation effect 
plays an important role in describing the charge carrier properties. 



IV. SUMMARY AND CONCLUSION 



For a model conjugated polymer chain initially holding a polaron defect, described using 
the combined SSH-EHM model extended with an additional part for the influence of an 
external electric field, we have studied the dynamics of polaron transport through this chain 
by virtue of simulating the backbone monomer displacements with classical MD and evolving 
the wavefunction for the vr-electrons with the adaptive TDDMRG. 

It is found that after a smooth turn-on of the external electric field the polaron is ac- 
celerated at first up to a stationary constant velocity as one entity consisting of both the 
charge and the lattice deformation. During the entire time evolution process the geomet- 
rical distortion curve and charge distribution shape for the polaron defect always stay well 
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coupled and show no dispersion under low external electric fields, implying that the polaron 
defect is an inherent feature and the fundamental charge carrier in conducting polymer. The 
dependence of the stationary velocity of a polaron Vp on the external electric field strength 
is also studied, and an ohmic region where Vp increases linearly with the field strength is 
found. Values beyond the speed of sound are achievable. This ohmic region extends ap- 
proximately from 3 mV/A to 9 mV/A for the case of ?7=2.0 eV and 1^=1.0 eV. Below 3 
mV/A the polaron velocity increases nonlinearly and in high external electric field strengths 
the polaron will become unstable and dissociate. Our calculated threshold field strength for 
polaron dissociation is around 10 mV/A for the case of f/=2.0 eV and V=1.Q eV. It is in 
good agreement with experimental results and more physically intuitive than previous SSH 
calculations which take no electron-electron interactions into account. 

The influence of electron-electron interactions (both the on-site Coulomb interactions U 
and the nearest-neighbor interactions V) on polaron transport are investigated in detail. 
In general, the increase of the on-site Coulomb interactions U makes the lattice tend to be 
occupied by one electron per site and accordingly suppress the polaron transport. Therefore, 
Vp decreases monotonically with U . Meanwhile, small V values are not beneflcial to the 
polaronic motion because they induce a more localized defect distribution, and due to the 
induced large charge polarization accompanied with large nearest-neighbor attractions large 
V values favor the polaron transport. When U and V are considered at the same time, Vp 
also decreases monotonically with the increasing U value and Vp decreases more rapidly for 
small U in large V cases than in small V cases. 
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